Analytic thermodynamics and thermometry of Gaudin-Yang Fermi gases 
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We study the thermodynamics of a one-dimensional attractive Fermi gas (the Gaudin-Yang model) 
with spin imbalance. The exact solution has been known from the thermodynamic Bethe ansatz for 
decades, but it involves an infinite number of coupled nonlinear integral equations whose physics is 
difficult to extract. Here the solution is analytically reduced to a simple, powerful set of four algebraic 
equations. The simplified equations become universal and exact in the experimental regime of strong 
interaction and relatively low temperature. Using the new formulation, we discuss the qualitative 
features of finite-temperature crossover and make quantitative predictions on the density profiles 
in traps. We propose a practical two-stage scheme to achieve accurate thermometry for a trapped 
spin-imbalanced Fermi gas. 



Alkali fcrmionic atoms at tens of nano-Kclvin are a re- 
markable addition to correlated quantum matter. Their 
interaction is unprecedentedly tunable from attractive to 
repulsive infinity. This allows, for instance, the observa- 
tion of universal properties of Fermi superfluids at uni- 
tarity. Recent developments [l], H} provide imbalanced 
populations of cold atoms in different hyperfine (spin) 
states, adding a new dimension in the low-temperature 
phase diagram In low dimensions, such attractive 
Fermi gases with spin imbalance give rise to a series of 
very interesting phases. For example, it was suggested 
that the quasi one-dimensional imbalanced Fermi gas, 
i.e., a weakly coupled array of one-dimensional (ID) at- 
tractive Fermi gases, gives a better chance to observe a 
long-sought crystalline superfluid, known as the Fuldc- 
Ferrell-Larkin-Ovchinnikov (FFLO) state [J, |6|. 

The realization of ID ultracold atomic gases 0, 0] pro- 
vides a new setting to experimentally test exactly solved 
models of interacting fermions. In particular, spin imbal- 
ance enables us to reach parameter regimes impossible 
in solid state (such as quantum wires) to study exotic 
quantum liquids outside the paradigm of a spin-charge 
separated Tomonaga-Luttinger liquid (TLL). Since their 
thermodynamics can be computed exactly, ID ultracold 
Fermi gases can also serve as a calibration reference sys- 
tem to measure the thermodynamics of strongly interact- 
ing Fermi gases in higher dimensions. 

Contrary to the impression one might have, ex- 
act "solvability" (integrability) does not guarantee that 
physical quantities of interest can be actually calculated. 
In principle, the thermodynamic Bethe ansatz (TBA) [9( 
gives exact results for thermodynamic properties. How- 
ever, it requires numerically solving an infinite number 
of coupled nonlinear integral equations to compare with 
experiments jT(|. On the other hand, quantum Monte 
Carlo calculations are limited to small system sizes 11 1. 

In this paper, we give a new formulation of the exact 



thermodynamics of a ID Fermi gas with contact attrac- 
tion, the Gaudin-Yang model 12j, Il3fl , which is now ac- 
cessible in cold atom experiments [3, @|. Based on the 
analytical analysis of the TBA equations, we derive a 
simple, complete set of algebraic equations which yield 
all thermodynamic quantities and show that they are 
asymptotically exact and universal in the physically in- 
teresting regime. Our approach can be applied to a wide 
range of Bethe ansatz integrable many-body systems. As 
an application, we propose a two-stage scheme to achieve 
accurate thermometry for trapped ID Fermi gases. 

Model and strong coupling limit. In the on-going ex- 
periments, for instance, at Rice University Q, a system 
of parallel ID gas "tubes" is prepared by loading alkali 
fermionic atoms (e.g., 6 Li in hyperfine states \F = 1/2, 
772 p = ±1/2}) in a square optical lattice. The trans- 
verse dynamics is suppressed by deep lattice potentials, 
and the motion is restricted along the tube axis x. The 
Fermi gas confined in each tube is then well described by 
the Gaudin-Yang model 

nana, 



2m dx 2 



1p<r(x) 



(1) 



The two hyperfine species with equal mass m are labeled 
with spin up and down respectively, g > is the contact 
attractive interaction. In general, the chemical poten- 
tials for the spin up and down fermions are different, 
say /Zf > Following convention, we define the chem- 
ical potential /i = (/ij + /i|)/2, the effective magnetic 
field ft, = (/ij — A*i)/2, the total density 77 = n-\ + Jij., 
the magnetization M = — rij_, and the polarization 
P = M/n. There are two length scales in this problem: 
the ID scattering length ai = 2ft 2 /(rag) characterizing 
the interaction strength g, and the inter-particle spac- 
ing 1/n. Their ratio defines the dimensionless interaction 
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strength 7 = 2/(na\). In the experiments, the gas is usu- 
ally dilute and strongly interacting, namely a\ <C 1/n, so 
7> 1. We will focus on this strong coupling limit where 
the binding energy Eb = h 2 /(ma 2 ) is much larger than 
the chemical potentials. 

The Gaudin-Yang model is exactly solvable by means 
of the Bethe ansatz 12 , 13| • Its zero temperatu re p hase 
diagram has been worked out theoretically [3, H, [3] • 
There are three phases: the fully paired (BCS) phase 
which is a quasi-condensate with zero polarization, the 
fully polarized (normal) phase with P = 1, and the par- 
tially polarized (ID FFLO) phase where < P < 1 and 
the pair correlation function oscillates in space |5|. For 
given /z, the FFLO phase is separated from the BCS 
phase and the normal phase by two quantum critical 
points at ft, = h c \ and /i C 2, respectively. An effective field 
theory of the ID FFLO phase, recently established by 
two of us, shows that it is a novel type of two-component 
TLL, each gapless normal mode being a hybridization of 
spin and charge [5(. 

Here, based on the exact Bethe ansatz solution, we 
demonstrate that, in the low-energy limit the ID FFLO 
is equivalent to two coupled gases of spinless fermions. 
Due to strong attraction, any minority (spin J,) fermion 
would like to pair up with a majority (spin |) fermion. 
Then, roughly speaking, the FFLO phase can be viewed 
as a mixture of tightly bound pairs (labeled by subscript 
b) and unpaired leftover fermions (labeled by subscript 
u) 0, M, 0. By a transmutation in statistics, the 
bosonic pair degree of freedom is equivalent to a gas of 
spinless fermions, analogous to the well known case of 
Tonks-Girardeau gas of hard-core bosons. The two gases 
are coupled by residue scattering, so the effective chemi- 
cal potential of each gas, (/J- u ), depends on the Fermi 
pressure of the other gas, p u (pb)- 

Exact thermodynamics. Our analysis begins with rec- 
ognizing the separation of energy scales in the strong 
coupling limit. At very low temperatures, T <C fib 
(Boltzmann's constant ks = 1), the thermodynamics of 
the ID FFLO phase is governed by the linearly dispersing 
phonon modes, i.e., the long wavelength density fluctua- 
tions of the two weakly coupled gases. Using the effective 
field theory and the standard conformal mapping [l7l.[T^. 
we find the low temperature thermodynamic potential 
(per unit length), 



f FT (T)=%-^-(- + -)T 2 
6ft v b v u 



(2) 



The system has two gapless excitations, which are the 
long wavelength density fluctuations (phonons) of the 
two gases. Their group velocities are the respective 
Fermi velocities, vj, = (1 — P)[l + (1 + 3P)/j]v wjA and 
v u = P[l + 4(1 - P)/~/]vp, with vf = hnir/m [HI]. At 
higher temperatures, T ~ fib or fi u , the particle and hole 
excitations are no longer restricted near the Fermi sur- 
face. Their dispersion undergoes a crossover from rela- 



tivistic (linear in k) to non-relativistic (fc 2 ). In particular, 
this happens in the vicinity of quantum critical points, 
where either fib or fi u becomes vanishingly small. The 
conformal invariance is now violated due to the lack of 
Lorentz invariance. Consequently, Eq. @ becomes inad- 
equate. Nevertheless, we shall show that the system is 
still equivalent to two weakly coupled gases despite the 
lack of conformal invariance. 

At even higher temperature, T ~ Sb, there is enough 
thermal energy to break pairs. Similarly, for T ~ 2h, a 
leftover fermion can flip its spin to be anti-parallel to the 
external field h due to thermal excitation. However, in 
this paper, we focus on the physically interesting regime 
of low temperatures and strong coupling, 

(3) 



T<C£ B ,2ft and 7>1. 



Our key observation is that the TBA equations can be 
greatly simplified in the regime ([3]) due to the suppres- 
sion of spin fluctuations. Mathematically, the infinite 
set of integral equations for spin rapidities is analytically 
tractable for strong attraction and high magnetic field. 
Their overall effect can be absorbed into a renormaliza- 
tion of fi u , which becomes exponentially small at low T 
so the feedback effect to the spin rapidity solution can 
be safely neglected. This leads to our central result — 
computing the thermodynamics at finite temperatures 
only requires solving the coupled algebraic equations: 



ax 



Vb = M + ~2 ~ -jPb 



aiPu, 



Mu = A* 
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Tln(2cosh-) 



aiPb 



T 3/2 Li3(-e 2rt 



2irh 2 



T^ 2 Us(-e^/ T ). 



(4) 

^p u sech 2 ^, (5) 
(6) 



(7) 



Here, Li s (x) is the standard polylogarithm function as a 
result of the Fermi-Dirac integral. As we shall focus on 
the regime ([3]), we can further neglect the exponentially 
small contributions from the spin fluctuation. Hence, 
Eq. ([5]) is replaced by 



axpb- 



(8) 



The thermodynamic potential (per unit length) & is 
then given by Sf = — p = —{pb+p u ), where p is the 
pressure. From Sf, it is straightforward to find the den- 
sity n — —d^/dfi, the magnetization M = —d^/dh, 
the entropy s — —d^/dT, and the compressibility k = 
dn/dfi. For these purposes, it is convenient to first 
take the corresponding derivative of Eqs. dD)-©, use 
<9Li s (— e x )/dx = Li s _!(— e x ), and then solve the re- 
sulting linear equations. For example, defining a — 
-(mT/4irh 2 ) 1 / 2 U 1/2 (~e 2 ^/ T ), (5 = -(mT , /2 7 rfc 2 ) 1/2 
Li 1 / 2 (— e M "/ T ) and c = 1 + a\a — 4a 2 a/3, we have 

n= (ia + /3-7aial3)/c, M = 0(1 — Za\a)/c. (9) 
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FIG. 1: (Color online) Left panel: Comparison of the tem- 
perature dependence of the pressure p obtained by numeri- 
cally solving the full TBA equations (circles) and simplified 
Eqs. (I!}-© (solid line). The entropy s(T) (green dashed 
line) is linear in T at low temperatures, agreeing well with 
sft (blue dotted line) as predicted by the low energy effec- 
tive field theory, /i = —0.95 and h = 1.1. The energy and 
length units are £s/2 and oi, respectively. Right panel: The 
schematic finite-temperature "phase diagram" in the regime 
Eq. ©. 



A detailed analysis of the TBA equations together with 
the derivation of Eqs. ©-© will be presented elsewhere. 

For given (//, h, T), equations ©-© can be solved by 
simple iteration by treating aipi, as being small compared 
to fit, and similarly a\p u as to fj, u . These equations free us 
from the requirement of solving the full TBA equations, 
which is a much harder task [10(. To test the accuracy 
of our new approach, we examined the pressure p = — §f . 
In Fig. [H we compare results from Eqs. ©-© (solid 
lines) and that obtained by numerically solving the full 
TBA equations (circles). We observe that indeed only 
at high temperature, T ~ £b/2, the deviation from the 
exact result becomes significant. Thus the validity of our 
new formulation is established. 

In the low temperature limit T <C ^ , /i„ we calculate 
the free energy using a Sommerfeld expansion based on 
Eqs. ©-©. Indeed, we find the thermodynamic poten- 
tial reproduces the field theory result [Eq. ©] and the 
corresponding entropy (per unit length) 



sft 



n 1 + 1 ) T 

3h v b v u 



(10) 



The entropy s(T) obtained from Eqs. I©-© is shown 
in Fig. [T] together with the low energy effective field 
theory result sft- s (T) is clearly linear in T at low 
temperatures and agrees well with sft- The deviation 
from the linear T dependence marks a universal crossover 
around T* from relativistic to non-relativistic dispersion 
for particle- and hole- like excitations. This leads to a 
minima in the magnetization as well as in the densities 
(not shown in Fig. © at roughly the same temperature 
scale T* . Analogous crossover phenomena were discussed 
in the mag netization curve for gapped spin-1 Heisenberg 
chains 11911 and experimentally observed in two-leg spin 
ladders [20j ]. 

The gross features of our findings are summarized in 



the schematic "phase diagram" in Fig. Q] There is no 
phase transition at finite temperature but only crossovers 
between different regimes. At low temperatures, the 
BCS, normal, and FFLO phase become the relativistic 
TLL of bound pairs (TLLf,), unpaired fermions (TLL U ), 
and two-component TLL (TLL2) respectively. The TLL 
description of the FFLO phase breaks down at crossover 
temperature T* where the dispersion of cither bound 
pairs or leftover fermions becomes non-relativistic. In 
particular, in the vicinity of the quantum critical points 
(h c i and h C 2), the system crosses over to the "quantum 
critical (QC) regime". The other crossover (dash-dot) 
lines roughly correspond to the excitation gap for un- 
paired fermions or bound pairs. 

Application: two-stage thermometry. Now we use these 
formulas to study the strongly interacting imbalanced 
Fermi gas inside a ID harmonic trap of frequency u>. Ac- 
cording to the local density approximation, the chemical 
potential varies in space as fx{x) = fio — ^muj 2 x 2 , with 
x = corresponding to the center of the trap, while 
the magnetic field h stays constant. Thus, distinct zero 
temperature phases are realized at different locations in 
the trap, giving rise to particular spatial structures in 
the in-situ density images. A central challenge in inter- 
preting the density distribution data from experiments is 
how to determine (T, fiQ,h), none of which seems to be 
accurately measurable so far. 

We use our results to propose a two-stage scheme to 
accurately determine the parameters T, fio, and h from 
the density profiles. Stage 1: estimate the values of these 
parameters using the density profiles at the phase bound- 
ary and at the edge of the gas cloud. Stage 2: using the 
estimations as guess inputs, obtain (T, Ho,h) to high ac- 
curacy by a direct three-parameter fit of the whole den- 
sity profile using Eq. ©. Previously, thermal tails at the 
outer wing of the density cloud have been used to extract 
the temperature of atomic gases in 3D traps based on 
the observation that the gas is essentially non-interacting 
there 2l|, |22j. Here, by contrast, the outer wing may still 
interact strongly (e.g., it might be a Tonks-Girardeau 
gas). Moreover, our method utilizes densities not only at 
the edge of the cloud, but also at the interface. We now 
illustrate how to accomplish the goals of Stage 1 for the 
cases of low and high total polarization. 

For low polarization, the trap consists of a partially 
polarized (FFLO) core and a fully paired (BCS) wing, 
which at strong coupling is a Tonks-Girardeau gas. Rep- 
resentative finite T density distributions are shown in 
the left panel of Fig. O At T = 0, the magnetiza- 
tion M drops to zero at the FFLO-BCS phase bound- 
ary /i = —h + n 3 7r 2 ai/24 [l6| . while 7if and vanish 
at fi — —sb/Z- Finite temperature leads to thermal 
tails for all three. The compressibility K is also shown 
in Fig. [21 which develops a peak near the edge of the 
cloud. At T = 0, k diverges as ~ 1/vV + £ b/2 as fi ap- 
proaches — £b/2. At finite temperature, this divergence 
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FIG. 2: The density nx, n\, the magnetization M — n-| - — n\, 
and the compressibility k = dn/dfj, versus the local chemi- 
cal potential /i in a ID trap. Left panel: low polarization, 
h = 0.9. Right panel: high polarization, h — 1.1. The tem- 
perature T = 0.02 (units are the same as Fig. [T}. k develops 
a peak at the thermal tail of the gas or at the phase boundary. 
Its peak value depends sensitively on temperature. 

is replaced by a peak which is increasingly shifted away 
from /i = — £b/2 and becomes less pronounced as the 
temperature is raised. In the BCS region, p u is essen- 
tially zero, and the equation of density [Eq. (JS])] simplifies 
to ri = 4a/ (1 + aid). An analytic formula for k follows 
with its peak value given by the simple expression 

, 1 2m 1.29 , , x 

^ = 0.712^—--— + ...]. (11) 

For strong interaction and low temperatures, the second 
and higher order terms can be neglected. This formula 
provides a powerful way to estimate T, namely the peak 
of compressibility is a sensitive thermometer. In fact, T, 
/Jo, and h can be determined from experimentally mea- 
sured n(x) and M(x), in the following steps. (A) Obtain 
the profile of compressibility from the density distribu- 
tion in the trap, k{x) = — (muj 2 x)~ x dn / dx . (B) Read off 
the position, aipk, and the magnitude of the compress- 
ibility peak, Kpk, from the plot k(x). (C) Calculate T 
from K p k, using Eq. (fTTj) . (D) Estimate the value of 
the chemical potential (Iq at the center of the trap from 
Vb(x P k) = Mo - muj 2 x 2 pk /2 + e B /2 ~ 0.554T. (E) Esti- 
mate h using h = — Mo + muj 2 x 2 M /2, where Xm is read 
off from the magnetization profile M(x) at the value of 
M(x M ) = 0.605 y^mT/2TrW. Note that the thermal tail 
of M is usually very broad, because [i u <C Mb is typically 
smaller than T for low polarization. 

For high polarization, a typical finite T density profile 
is shown in the right panel of Fig. [2] At zero temper- 
ature, the trap consists of a partially polarized (FFLO) 
core and fully polarized normal (N) wing, with n± de- 
caying to zero at the FFLO-N phase boundary. Thermal 
fluctuations smear the T — phase boundary, but the 
compressibility develops a pronounced peak at 
due to the rapid decay of ri| at low temperatures. In 
the normal region at the outer wing, the fully polarized 
gas is free and we found that the compressibility has the 



simple form k = -(m/27rfi 2 T) 1 / 2 Li_ 1/2 (-e' i " /T ). In this 
region, the thermal tail of n-f leads to a second peak in 
the compressibility at x = x V 2- The two-peak structure 
of the compressibility is clearly seen in Fig. O The first 
stage estimation of (T,fj,o,h) in this case involves three 
steps. (A') Measure the value of compressibility, k P 2, 
at the second (outer) peak, and obtain the temperature 
from the relation 

« p2 = 0.12^ — . (12) 

(B') Estimate no from the first (inner) peak location 
x p i, using fib(xpi) — fio — mu) 2 x 2 1 /2 + eb/2 ~ 0.554T. 
(C) Estimate h from the second peak location x — x P 2 
using fJ, u (x P 2) = Mo _ muj2x l2 + h= 1.1 IT. 
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